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1. Introduction 

Massive stars, losing energy to radiation and photons, evolve until an iron core is 
formed (for a recent review see Ref. ^). This core, which has a very low entropy 
per baryon, is supported by the electron degeneracy pressure. As a consequence such a 
core is dynamically unstable and collapses until the matter is mostly neutronized and 
supernuclear densities are reached. Only when the density significantly exceeds the 
nuclear density, the pressure becomes sufficiently repulsive to stop the collapse. In the 
current paradigm the innermost shell of matter reaches to these densities first, rebounds 
and sends a pressure wave through the rest of the core. Such waves, produced by the 
subsequent shells and traveling faster than the infalling matter, collect near the sonic 
point. As that point reaches nuclear density a shock wave breaks out. The subsequent 
evolution of this bounce shock is not yet well-understood and is subject to much study 
13 El- Current models fail to explode. 

Even though it has not yet been demonstrated that explosion is an outcome of the 
core-collapse, it is well-established that the newly-formed hot proto-neutron star cools 
by neutrino emission. (It was those neutrinos that were observed in Supernova 1987A). 
Almost all (99%) of the gravitational binding energy of the neutron star 
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is radiated away in neutrinos of all flavors. Altogether a star with mass ~ 8Mq will emit 
~ 10^^ neutrinos. Thus, after the explosion ejects the material from the outer layers, a 
"neutrino-driven" wind may blow the medium above the neutron star, heating it to an 
entropy per baryon of several hundreds in units of Boltzmann's constant. For such high 
entropies nuclear statistical equilibrium is not established for nuclei heavier than alpha 
particles. 

The rapid neutron-capture process (r-process) is responsible for the formation of a 
number of nuclei heavier than iron. (For a recent review see Ref. Q). The astrophysical 
site of the r-process nucleosynthesis is not yet identified. For r-process nucleosynthesis 
to successfully take place a large number of the neutrons are required to interact in 
a relatively short time, indicating that r-process sites are associated with explosive 
phenomena. Indeed the seminal Burbidge, Burbidge, Fowler, and Hoyle paper suggested 
the neutron-rich ejecta outside the core in a type II supernova as a possible site of the r- 
process More recent work pointed to the neutrino-driven wind in the supernovae as a 
possible site [El 13 IE] • Meteoric data and observations of metal-poor stars indicate that r- 
process nuclei may be coming from diverse sources jH] . Binary neutron star systems were 
also proposed as a site of the r-process (see e.g. Ref [ID])- In outflow models r-process 
nucleosynthesis results from the freeze-out from nuclear statistical equilibrium. The 
outcome of the freeze-out process is determined by the neutron-to-seed ratio. This ratio 
in the post-core-bounce supernova environment is controlled by the intense neutrino 
flux radiating from the neutron star. 

Neutrino interactions play a crucial role in core-collapse supernovae. (For a brief 
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summary see Ref . ^T] ) . Neutrino heating is one of the possible mechanisms for reheating 
the stalled shock ^21- The neutrino fluxes control the proton-to-neutron ratio in the 
high-entropy hot bubble. As we describe in the next section there is a hierarchy of 
energies for different neutrino flavors. Hence swapping active neutrinos via neutrino 
oscillations changes the n/p ratio and may alter r-process nucleosynthesis conditions 
[T^ . Neutrino oscillations in a core-collapse supernova differ from the matter-enhanced 
neutrino oscillations in the Sun as in the former there are additional effects coming from 
both neutrino- neutrino scattering fT^ ITH] and antineutrino flavor transformations ^Hl • 
We present a simple description of the core-collapse supernovae in the next Section. 
In this section after summarizing properties of the neutrino-driven wind, we discuss 
neutrino fluxes, luminosities, reaction rates, and the equilibrium electron fraction. A 
brief description of neutrino mixing and neutrino interactions relevant to core-collapse 
supernovae is given in Section 3. In Section 4 we calculate values of the electron fraction 
under several evolution scenarios that may impact r-process nucleosynthesis. 



2. A Simple Description of Core Collapse Supernovae 

2.1. Neutrino-Driven Wind 

A careful treatment of the neutrino-driven wind in post-core bounce supernova 
environment was given in Ref. Here we present a heuristic description following 

Ref. PH]. One can assume that at sufficiently large radius above the heating regime 
there is hydrostatic equilibrium |19j : 

dP ^ GMnsP ^2) 
dr ' 

where P is the hydrostatic pressure, G is Newton's constant, M^s is the mass of the hot 
proto-neutron star, and p is the matter density. Using the thermodynamic relation for 
the entropy at constant chemical potential, /x, 

5'total = ( 7777 ) ) (3) 



6T J 



and integrating Eq. (0) we can write entropy per baryon, S, as 

TS = (4) 

r 

where is the average mass of one baryon, which we take to be the nucleon mass. In 
the region above the neutron star the material is radiation dominated and the entropy 
per baryon can be written in the relativistic limit as 

k 45 pB\hcJ ' 
where the statistical weight factor is given by 

9s = 3b + l 3f- (6) 

bosons fermions 
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Assuming a constant entropy per baryon, Eqs. (jH) and (jH)) give the baryon density, ps, 
in units of 10^ g cm~'^ as 

where S'loo is the entropy per baryon in units of 100 times Boltzmann's constant, Vj is 
the distance from the center in units of 10^ cm, and we assumed that Mns = 1.4Mq. 
Defining Tg to be temperature in units of 10^ K, Eq. takes the form 

Tg^ioo ~ • (8) 

In Figure d we present matter density and temperature profiles based on heuristic 
description given in this section. Several values of S'loo can be used to describe stages 
in the evolution of supernovae. Smaller entropies per baryon, Siqq < 0.5, provide a 
better description of shock re- heating epoch, while larger values, Sioo ^ 1, describe late 
times in supernova evolution, namely, the neutrino- driven wind epoch. Higher entropy 
corresponds to less ordered configurations with smaller baryon densities. In Figure 1 the 
statistical weight factor, gg, is taken to be 11/2 in the calculation of matter density since 
temperature and entropy per baryon are Tg > 4 and S'loo ^ 1-5, respectively. Under 
these conditions both photons are electron-positron pairs are present in the plasma. 
When temperature drops, Tg < 4, only photons present in the medium and statistical 
weight factor has to be taken gg ^ 2. 



2.2. Neutrino Fluxes and Luminosities 



We adopt the prescription in Ref. J4] for the neutrino fluxes. If we take the density 

of particles to be p, the number of particles that go through an expanding surface of 

radius r per unit time is 

dN dV ^ r, dr , , 

_ = = i,r-,- (9) 

where is the total number of particles and V = Anr^. Assuming a Fermi-Dirac 
distribution function for the neutrino number densities, ignoring neutrino mass in 
comparison to its energy, and taking the particle velocity to be c, Eq. gives for 
the flux of neutrinos emitted from the neutrinosphere as 

d^^ 



^'"^^^ dn,, (10) 



87r3 (nc)3l + exp[(E,-/i,)/A;T, 
where /i,^ is the neutrino chemical potential. Choosing the z-direction as the vector that 
connects the point where the flux is to be calculated (with radial position r) to the 
center of the neutron star, one sees that the azimuthal symmetry still holds, but the 
polar angle is bounded by the finite size of the neutrinosphere. Thus 

cil]^ = 27r(l-cos^o), (U) 
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Figure 1. Solid (5ioo = 0.5), dashed (5*100 — 1), and dotted (5ioo=l-5) lines 
corresponds to matter density profiles (upper panel), and temperature profiles (lower 
panel) based on heuristic description. The thick band in the upper panel is matter 
density profile from numerical supernova models for tpg w 4 s (taken from ref |13|). 



where 



i?2 i?2 

^0^^0 = ^/1-3^-1-^ (12) 



r 



with being the radius of the neutrinosphere. Within this approximation one then 
obtains the differential neutrino flux as 

d^u _ I c Rl El ^^g^ 



dE^ 87r2 {hcf r2 1 + exp[(E^ - n^)/kT^] 



Neutrino Mixing and Nucleosynthesis in Core-Collapse Supernovae 



6 



Using similar reasoning one can write an expression for the neutrino luminosity. 
Replacing N in Eq. Q by the total energy and the matter density p by the energy 
density one can write down 

2 1 /" Eyd^Y*u 

= "^2^ J l + exp[{E,-p,)/kny ^^^^ 

Again ignoring the neutrino mass as compared to its energy and using Eqs. (fTT|) and 
(fT^ to do the angular integration we obtain 

where rj = p^/kTy and F^{ri) is the relativistic Fermi integral 

/•oo 3 

In our calculations we take 77 = and use the value -^3(0) = 77r'^/120. Often it is 
convenient to express the neutrinosphere radius in terms of the neutrino luminosity. 
Rewriting R,^ in terms of Ly and inserting it into Eq. (jl3j) we get 



( El 



dE, 47rr2 {kT,YF^{r]) \l + exp[(E, - p^)/kT, 
2.3. Reaction Rates 



(17) 



The dominant reactions that control the n/p ratio is the capture reactions on free 
nucleons 

z/e + n^p + e~, (18) 

and 

z/e + P ^ n + e+. (19) 
We take the cross sections for the forward reactions to be ^H] 

a.XE..) - 9.6 X 10-^^ ( ^tv"' )''"^'' 

and 

a.,(E,J ^ 9.6 X 10-^^ ( ^'.v" )'^"^'' ^2^^ 

where A„p 1.293 MeV is the neutron proton mass difference. For simplicity we 
ignored weak magnetism and recoil corrections, which may be important [20] • These 
corrections cancel for the former cross section, but add for the latter one. The rates of 
these reactions can be written as 

To calculate the neutrino capture rates on nuclei one needs to include all possible 
transitions from the parent to the daughter nucleus, including not only the allowed 



A = / a{E),^dE,. (22) 



Neutrino Mixing and Nucleosynthesis in Core-Collapse Supernovae 



7 



transitions, but also transitions to the isobaric analog states, Gamow- Teller resonance 
states, transitions into continuum, as well as forbidden transitions. Aspects of such 
calculations are discussed in Refs. [22] and [2n! (see also ^|). Careful input of such 
reaction rates in supernova simulations is especially crucial to assess the possibility of 
core-collapse supernovae as a site of the r-process nucleosynthesis [HI [2^] • 

Because of their charged-current interaction electron neutrinos may play a role in 
reheating the stalled shock as well as regulating the neutron-to-proton ratio. In contrast, 
since the energies of the muon and tau neutrinos and antineutrinos produced are too 
low to produce charged leptons, these neutrinos interact only with the neutral-current 
interactions. Recently significant attention was directed towards understanding the z/^ 
and spectra formation. Neutrinos remain in local thermal equilibrium as long as they 
can participate in reactions that allow exchange of energy and neutrino pair creation or 
annihilation. It turns out that the neutrino bremsstrahlung process 



is more effective than the annihilation process u + u ^ e'^ + e~ at equilibrating 
neutrino number density [2^1 • (However the neutrino- neutrino annihilation process 
Ue + Ue ^ i^/i + i^fi is one of the primary sources of muon and tau neutrinos j28j). 
The impact of the neutrino bremsstrahlung process on equilibrating the energy spectra 
seems to be comparable to that of 



which dominates the neutrino spectra formation. Finally recoil corrections to the uN 
interactions are very important in the formation of the z/^ and Ur spectra as they permit 
energy exchange |29j . 

2.4- Electron Fraction 

The electron fraction, Yg, is the net number of electrons (number of electrons minus the 
number of positrons) per baryon: 



where Ug-, ne+, and ub are number densities of electrons, positrons, and baryons, 
respectively. Introducing A^^, number of species of kind j per unit volume, and Aj, 
atomic weight of the j-th species, one can write down expressions for the mass fraction. 



N + N^N + N + u + iy 



(23) 



The most effective process to exchange energy is fI7\ 
u + N + N^N + N + u, 



(24) 



(25) 




(26) 




(27) 



and the number abundance relative to baryons, Yj 
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The electron fraction defined in Eq. (|26|) can then be rewritten as 

i i ^ ^ 

= Xp + Jx„ + V (^] Xh, (29) 



h ^ 

where Zi is the charge of the species of kind i, and the mass fractions of protons, Xp, 
alpha particles, X^, and heavier nuclei ("metals"), X^, are explicitly indicated. 
The rate of change of the number of protons can be expressed as 

^ = -(Ap, + K-)N, + (A,, + Ae+)X„, (30) 

where X^^ and Ag- are the rates of the forward and backward reactions in Eq. ()18|) and 
Ap^ and Ae+ are the rates of the forward and backward reactions in Eq. (fTTHl . Since the 
quantity ^ • NiAi does not change with neutrino interactions, one can rewrite Eq. (jHUjl 
in terms of mass fractions 
dX 

= -(Ap, + Ae-)Xp + (A,, + Ae+)X„. (31) 

In the hot bubble the rates are usually expressed in terms of the radial velocity field, 
f (r), above the neutron star, i.e. dY/dt = v{r)[dY/dr]. A careful study of the influence 
of nuclear composition on Ye in the post-core bounce supernova environment is given in 
Ref. pUj to which we refer the reader for further details. 
If no heavy nuclei are present we can write 

Because of the very large binding energy, the rate of alpha particle interactions with 
neutrinos is nearly zero and we can write dV^/dt = dXp/dt. Using the constraint 
Xp + X„ + Xq, = 1 and Eq. (jH^ . Eq. (pTTjl can be rewritten as 

rfFe 1 



An — (Ap + A„)Ye + -(Ap — A„)X(^, (33) 

where we introduced the total proton loss rate Ap = Ap^ + Ae- and the total neutron 
loss rate A„ = X^^ + Ae+. It has been shown that when the rates of these processes are 
rapid as compared to the outflow rate a "weak chemical equilibrium" is established |13j . 
The weak freeze-out radius is defined to be where the neutron-to-proton conversion rate 
is less than the outflow rate of the material. If the plasma reaches a weak equilibrium 
stage then Yg is no longer changing: dY^/dt = 0. From Eq. ()33p one can write the 
equilibrium value of the electron fraction 

Ap + A„ 2 Ap + Xn 

Different flavors of neutrinos decouple at different radii. Since z/^ and (and 
their antiparticles) interact with the ordinary matter only with the neutral current 
interactions, they decouple deeper in the core and have a large average energy. Electron 
neutrinos and antineutrinos have additional charged-current interactions with neutrons 
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and protons respectively. Since in the supernova environment there are more neutrons, 
electron antineutrinos decouple after z/^'s and z/^'s, but before electron neutrinos. 
Consequently one has a hierarchy of average neutrino energies: 

{E^J < {E,J < (35) 

where stands for any combination of i/^'s and i/^-'s. However a more complete 
description of the microphysics suggests that this hierarchy of average energies may 
not be very pronounced fl^ |2H1 1^ ■ This microphysics is dominated by the inelastic 
neutrino-nucleon interactions discussed in Section 2.3. 

In Figure |2l we present initial differential neutrino fluxes. In the upper panel the 
typical post-bounce neutrino energies are taken as the representative values of {E^J = 10 
MeV, {EpJ = 15 MeV, and {E^^^pJ = 24 MeV. (Average neutrino temperature for each 
flavor can be calculated using the relation: T^, = {E^) /3.151). Here neutrino and 
antineutrino luminosities are taken to be equal for all flavors. We express the neutrino 
luminosities in more convenient units of 10^^ ergs s^^. Equal luminosities of L^^ = 1 
is typically a good approximation for the neutrino-driven wind epoch. However, at 
earlier epochs, z/g and z/g luminosities can be as large as L^^ = 10. Except, through the 
neutronization burst, for few milliseconds luminosity can reach values of Lf} = 100, 
an order of magnitude larger than z/g during the same period. In the lower panel of 
Figure 2 we examine a less-pronounced hierarchy of neutrino energies. Here we adopt 
the representative values of {E^J = 13 MeV, (EpJ = 15 MeV, and {E^,,pJ = 16 MeV. 
Here the luminosities of z/^, and are taken to be the half of the value of the equal 
luminosities of z/g and z/g. 

In our calculations in Section 4, we adapt these initial distributions of z/g, z/g, 
and Ux at neutrinosphere, i?j,=10 km, with the indicated luminosities and follow the 
evolution of differential neutrino fluxes (number of neutrinos per unit energy per unit 
volume). 



2.5. Alpha Effect 

At high temperatures alpha particles are absent and the second term in Eq. ()34|1 can be 
dropped. In the region just below where the alpha particles are formed approximately 
one second after the bounce, the temperature is less then ~ 1 MeV. Here both the 
electron and positron capture rates are very small and Y^. can be approximated as 

As the alpha particle mass fraction increases (when Tg drops below 8) free nucleons 
get bound in alphas and, because of the large binding energy of the alpha particle, cease 
interacting with neutrinos. This phenomenon is called "alpha effect" [22] • Using Eq. 
(fHU)) one can rewrite Eq. (jH^ as 

Ye = Y}'^ +(1- Y}'A X„. (37) 
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E (MeV) 

Figure 2. Initial differential neutrino fluxes in arbitrary units. Solid, dashed, dotted, 
and thick lines correspond to the distributions of z^e, i^e, and v^. a) Fluxes with 
[E^^) = 10 MeV, {E^^) = 15 MeV, and (-E^^.pJ = 24 MeV. Neutrino and antineutrino 
luminosities are taken to be equal for all flavors. Distributions of Vx,Vx have much 
longer tails up to 100 MeV not shown in the flgure. b) Fluxes with {E^^) = 13 MeV, 
(E'p^) = 15 MeV, and (-Ei/^.p^) = 16 MeV. and luminosities are taken to be equal. 
Luminosities of other flavors are taken to be half of the v^. and Pe luminosities. 

Hence if the initial electron fraction is small {y!P^ < 1/2) the alpha effect increases 
the value of Y^. Since higher implies fewer free neutrons, the alpha-effect negatively 
impacts r-process nucleosynthesis 

Neutrino oscillations, since they can swap energies of different flavors, can effect the 
energy-dependent rates in Eqs. (jHHjl and (jSZ|), changing the electron fraction. Indeed, 
transformations between active flavors heat up i^e's and increase driving the electron 
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fraction to rather large values. Consequently a very large mixing between active neutrino 
flavors would have prohibited r-process nucleosynthesis in a core-collapse supernova 

[nun!. 

Actually electron neutrinos radiated from the proto- neutron stars are just too 
energetic to prevent the alpha effect in most cases. One possibility to reduce Xi^^ is 
to convert active electron neutrinos into sterile ones which do not contribute to this 
rate. This possibility is explored in references [TH], and [HH] . 

3. Neutrino Mixing 

Neutrino interactions in matter is a rich subject (for a brief review see Ref. |21])- While 
neutrinos are produced through weak interactions in flavor eigenstates, they propagate 
in mass eigenstates. Mixing angles correspond to rotations describing unitary connection 
between two bases: 









( 


\ 














J 









In our discussion, we follow the notation in Ref. [23 and denote the neutrino mixing 
matrix by Uai where a denotes the flavor index and i denotes the mass index: 

/I \ / Ci3 Sl,\ I Ci2 5i2 0\ 

= ^23 ^23 1 -Si2 Ci2 . (39) 

\0 -^23 ^23/ \-5i3 On) V 1/ 

In Eq. ()39j) C13, etc. is the short-hand notation for cos^^ia, etc. The notation Sl^, was 
used to indicate (sin 6'i3)e*'^ where is a CP- violating phase. We will ignore this phase 
in our discussion. 

We have compelling evidence supporting non-zero neutrino masses and mixings. 
Two-flavor solar neutrino solution corresponding to 6'i2 ~ 7r/6 and 5m^2 ~ 8 x 10~^ 
eV^ was identified using the recent Sudbury Neutrino Observatory (SNO) results (cf. 
Refs. inn], |SZ1)- Measurement of anti- neutrinos from nuclear power reactors in Japan 
by the KamLAND experiment confirmed this solution and improved the limits on the 
solar mass square difference, ^ml2, significantly |SH]- While it is known that the solar 
neutrino mixing angle is large, but not maximal, the atmospheric mixing angle is large 
and could be maximal, ^23 ~ 7r/4, as shown by the Super-Kamiokande experiment |39j . 
The latter angle is consistent with the KEK-to-Kamioka oscillation experiment, K2K 
|4()j . The corresponding atmospheric mass square difference is 5m23 ~ 3 x 10~^ eV^. 

The size of the last mixing angle, ^13 is currently best limited by the combined 
CHOOZ IH] and Palo Verde 1^ reactor experiments and SK atmospheric data. This 
is due to the null results from the reactor i/g disappearance over the Sm^^ distance 
scale. The upper bound on this angle from KamLAND and the solar neutrino data gets 
stronger (especially in the region with small atmospheric mass square difference where 
CHOOZ reactor bound is relatively weak), and even dominates, as this data get refined 
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Both measurements of the width of the Z boson and oscillation interpretation 
of the neutrino data, with the notable exception of the LSND signal 01], favor three 
generations of light active neutrino species. If LSND is confirmed, the most likely 
explanation could be the existence of a fourth neutrino (a relatively heavier sterile 
neutrino) since agreement between KamLAND and combined solar experiments already 
disfavors the interpretation of the LSND anomaly with CP violation. 

In the region above the supernova core density is still high, but steeply decreases. 
Matter-enhanced oscillations mediated by the solar mass square difference are impossible 
in the region close to the core which is suitable for the r-process. However, at late 
neutrino driven epoch baryon density could be low enough to allow resonances through 
Smf^, which is comparable to the atmospheric mass square difference. At such late 
times neutrino flux is expected to be low. We examine prospects of r-process at such 
environments in the next section. 

We describe neutrino mixing within the density matrix formalism E^l HH} 14 7j . 
We assume that the electron neutrino mixes with a linear combination of /x and r 
neutrinos. The neutrino oscillations are mediated by the matter mixing angle for 
transformation between 1st and 3nd mass eigenstates, ^13, and the corresponding 
atmospheric mass square difference, Sml^. The mixing between fi and r neutrino flavors 
does not have any significant effect on our results as long as total luminosities and 
corresponding average energies are equal and their mixing is maximal. In this limit, 
mixing between 2nd and 3nd mass eigenstates, can be rotated away and effectively 
electron neutrinos oscillates into some linear combination of /i and r neutrinos j^. 

In the rest of the paper 6 and Sm"^ refer to ^13 and ^m^g. This assumption simplifies 
the discussion and allows us to write the two-flavor density matrices as 

p= I ) =i(Po + P-a) (40) 

\ Pxe Pxx I 

and 

= ^(Po + P-a), (41) 

where we introduced the polarization vectors for neutrinos and anti-neutrinos, Pp and 
Pp. The diagonal elements in these expressions are initially given by the expression in 
Eq. p7|) . Non-diagonal elements are initially zero but may become non-zero during the 
neutrino evolution. 

Equations governing the evolution of neutrinos and antineutrinos can be cast into 
the forms jlH] 




9,P„= <^+A„ + V2G1 



p 

and 

9,Pp= <'-Ap + V2Gp 



/ rfq(l-^)(Pp-P,) 



iVeZ+ / rfq(l-^)(Pp-P,, 



xPp,(42) 



xPp,(43) 
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where 

Ap = ^ (sin 2^x - cos 2^z) , (44) 

Integrating Eqs. (jl^ and ()43|) exactly in the supernova environment is, at the moment, 
an unsolved problem. Indeed the exact solutions of these coupled, non-linear differential 
equations are expected to be very complicated. Instead we adopt the approximation 
proposed in Ref. [H] and also adopted in Ref. jlOj. In this approximation one uses 
flux-averaged values to obtain 

d,Pp = (^+Ap + v^GpiVeZ + V2GFF{r){3 - J)) x Pp (45) 

and 

drPp = (-Ap + v^GpA^eZ + V2GFF{r){J - J)) x Pp (46) 

In these equations J is the polarization integrated over all momentum modes, and 
F{r) = i[l — (1 — Rl/r'^Y^'^] is the geometrical factor introduced earlier (cf. Eqs. (|TT|l 
and (PI). 



4. Results and Discussion 



In our calculations, we concentrate on the late neutrino-driven wind epoch which is 
expected to have larger entropy. We adopt 5*100 = 1-5 and the gs = 2, since in this 
epoch the temperature is too low at later times to have electron-positron pairs to be 
adequately represented. 

In Figure 01 we present differential neutrino fluxes at several stages of the evolution 
in order to provide a better understanding of this mechanism. The mixing parameters 
are chosen as ~ tt/IO and ~ 3 x 10~^ eV^. Solid, dashed, dotted, and thick lines 
correspond to the distributions of z/g, z/g, z/a;,and Ux- In each panel, dependence is 
removed. For this figure the initial distributions are taken as the values given in Figure 
2a). Columns are for L^^ = 0.001, 0.1, and 50 from left to the right corresponding 
to very weak, moderate, and very strong neutrino self interaction contributions to the 
evolution. Rows show neutrino flux at r = 75 km, 100 km, and 150 km, exhibiting the 
evolution of neutrino distributions along the neutrino path. These should be compared 
to the initial (at ~ 10 km) neutrino distributions given in Figure Et')- Electron fractions 
corresponding to these luminosities are given in the left column of Figure |31 

• First column of Figure El L^^ = 0.001: Initially differential neutrino fluxes 
for and i>x are the same at ~ 20 MeV as seen in Figure Efe-)- At energies lower than 
20 MeV, differential flux of is higher than that of I'x whereas at energies higher than 
20 MeV the numbers reverse. The transformation starts at the low energy tail of the 
distributions. As we gradually move to lower densities resonance moves up to higher 
energies. Until we reach the resonance at 20 MeV the luminosity of z/g's decreases. 
At r = 100 km electron neutrinos below 20 MeV are mostly swapped with u^s while 
above 20 MeV there is yet no significant transformation. This point corresponds to the 
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Figure 3. Evolution of differential neutrino fluxes in arbitrary units. Solid, dashed, 
dotted, and thick lines correspond to the distributions of Ve, i^e f/j, and D^^. In each 
panel, dependence is removed. Columns are for L^^ — 0.001, 0.1, and 50 from 
left to the right corresponding to very weak, moderate, and very strong neutrino self 
interaction contributions to the evolution. Rows are calculated at r = 75 km, 100 km, 
and 150 km showing the evolution of distributions along neutrino path. See text for 
details. 



minimum of A of Eq. ()22|) for neutrinos (after the 1/r^-dependence is taken out). Since 
antineutrinos are not transformed A is constant for them. As a resuhs, r ~ 100 km 
corresponds to a dip in Y^. (sohd hne in figure |3]( a)). Since above 20 MeV the initial 
differential flux of u^s is smaller than that of u^^s, transformations after r ~ 100 km 
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r (km) r (km) 

Figure 4. Initial neutrino fluxes and luminosities are taken to be those in Figure 
2(a) in the left column and those in Figure 2(b) in the right column. Sioo is taken to 
be 1.5. In panels (a) and (b), solid, dashed and dotted lines correspond to the equal 
luminosities of L^^ = 0.001, 0.1, and 50 for all flavors, indicating very weak, moderate, 
and very-strong neutrino self interaction contributions to the evolution. In panels (d), 
and (e) i^e and De luminosities are taken to be L^^ = 0.002, 0.2, and 200 (solid, dashed 
and dotted lines, respectively). In (d) and (e) the luminosities of other flavors are taken 
to be L^^ = 0.001, 0.1, and 100 (sohd, dashed and dotted hues, respectively), (a, d) 
Equilibrium electron fraction as a function of the distance from the core with mixing 
parameters 6*13 ^ tt/IO and drnl^ ~ 3 x 10^"^ eV^. (b, e) Equilibrium as a function 
of the distance from the core with mixing parameters 6113 ~ 7r/20 and Smf^ ~ 3 x 10""^ 
eV^ . (c) Same as (b) but the impact of alpha particle formation is included according 
to the Eq. H34I) . Yg is shown when Xa~0, 0.3 and 0.5 (thin, medium and thick lines) 
for L^^=0.001 and 50 (solid and dotted sets of lines), (f) Same as (e) but the impact 
of alpha particle formation is included as in (c). Yg is shown when Xa=0, 0.3 and 0.5 
(thin, medium and thick lines) for L^^ =0.002 and 200 (solid and dotted sets of lines). 



would increase the number of high energy u^s. As neutrinos travel to regions further 
away (r= 150 km, where baryon density is much lower), swapping of high energy tail 
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of the z/g and distributions is completed, and approaches to the asymptotic fixed 
value. 

• Second column of Figure El L^^ = 0.1: This column describes the 
same evolution of neutrinos except that neutrino self-interaction effects play a more 
pronounced role. The chosen value of the luminosity could be representative of such 
late times in neutrino-driven wind epoch. Because of the small contributions of the self 
interactions, resonance region is relatively wider. No dip is observed in Y^. (dashed line 
in figureE](a)) at r = 100 km because transformation of low and high energy ends of the 
distributions happen almost simultaneously. There is also small transformation between 
z/g's and Ux^s. 

• Third column of Figure |^ L^^ = 50: This is an extreme case to illustrate 
the limit at which neutrino self interactions dominate. Swapping of both neutrinos 
and antineutrinos occur and Ye reaches its equilibrium value rapidly. After the 
transformation, electron neutrinos and antineutrinos assume similar luminosities and 
distributions since they are swapped with u^s and z/^^'s. Note that this equilibrium 
value is again different from 0.5 because of the threshold effects (due to the neutron- 
proton mass difference) in the reaction cross sections (cf. Eqs. fl20|l and ()2Hl ). 

In the left column of Figure 4 initial neutrino fluxes and luminosities are taken to 
be those in Figure 2(a) and in the right column initial neutrino fluxes and luminosities 
are taken to be those in Figure 2(b). In Figure |^a) we present the equilibrium electron 
fraction, Y^., as a function of the distance from the core for the three different cases of 
neutrino flux, L^^ = 0.001, 0.1, and 50, each corresponding to one of the columns in 
Figure 01 In Figure Ufa) we use the same neutrino parameters as in Figure El One 
could argue that the relatively large mixing angle of vr/lO is already disfavored by 
CHOOZ. To explore the implications of a smaller mixing angle, in Figure j^b) we show 
electron fractions calculated with the more realistic mixing parameters 6 ~ 7r/20 and 
Sm"^ ~ 3 X 10^^ eV^. The dip in Ye plot is sharper since resonance region will be much 
narrower with the smaller mixing angle. 

In Figures m a) and b), to calculate the equilibrium electron fraction, we ignored 
possible effects of the alpha particles and used Eq. (jH^ to evaluate Ye. InjU^c), we 
explore possible effects of alpha particles by using Eq. pljl to calculate Ye for three 
different values of X^. The "alpha effect" is manifest: as Xa gets larger, it pulls the 
value of the electron fraction closer to 0.5. 

In the right column of Figure 4 initial neutrino fluxes and luminosities are taken 
as the generic "almost equal" energies and different luminosities case shown in Figure 
2(b). It is immediately obvious that the situation is markedly different in this case. Even 
though the initial z/g and z/g spectra are about the same (cf. Figure 2(b)), Ye initially 
is greater than 0.5 because of the effect of the neutron-proton mass difference in Eqs. 
(1201 and (j2H). When the effects of the neutrino self- interaction terms are minimal (the 
low-luminosity case indicated by the solid lines in (d) and (e)), after z/g's and z/^-'s swap 
z/g luminosity decreases signiflcantly causing a big drop in the value of Yg. When the 
effects of the self-interaction terms are dominant (the high luminosity cases indicated 
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by the dotted lines in (d) and (e)), all flavors swap and the electron fraction eventually 
reaches to the asymptotic values. The impact of alpha formation, when the hierarchy 
of average neutrino energies is less pronounced, is presented in Figure 4(f). 

5. Conclusions 

In this article we investigated conditions for the r-process nucleosynthesis at late 
neutrino-driven wind epoch in a core-collapse supernova. We considered the region 
where the matter-enhanced neutrino transformation is driven by Sm^^, which we took to 
be comparable to the mass difference observed in the atmospheric neutrino oscillations. 

We found that, when initial luminosities are taken to be equal with a pronounced 
hierarchy of neutrino energies, the asymptotic value (at large distances from the core) 
of the electron fraction always exceeds 0.5, hindering the r-process nucleosynthesis. In 
this case we found that neutrino self-interactions decrease the electron fraction. In 
general these self-interaction terms, unhke the MSW effect, tend to transform both 
neutrinos and antineutrinos. Hence one expects that when the self interaction terms are 
dominant (e.g when the neutrino luminosities are very large) the electron fraction would 
reach to the value of 0.5. However we showed that, because of the threshold effects 
in the neutrino interactions, electron fraction exceeds 0.5 even when the background 
neutrinos arc numerous. Clearly these conditions are not favorable for the r-process 
nucleosynthesis. We also found that, when background effects are small, with the 
parameters we adopted for baryon density and neutrino mixing there exists a region 
around 100 km in which electron fraction is smaller than 0.5. Even in this region alpha 
effect pulls the value of the electron fraction closer to 0.5. Although the conditions 
are favorable, the r-process in this region could not contribute reasonable quantities of 
elements since the baryon density is rather low. 

In contrast, when the initial luminosities of the u^s and i^e's are taken to be 
twice of the other flavors with a much less pronounced neutrino hierarchy, we found 
that the asymptotic electron fraction is less than 0.5 when neutrino self-interactions 
are negligible. However, once again the baryon density is low and r-process in this 
asymptotic region is likely to produce insigniflcant quantities. When neutrino self- 
interactions are dominant the asymptotic electron fraction shows the same behavior 
as the first generic case (with a pronounced hierarchy of neutrino energies and equal 
luminosities) we considered. 

It is not surprising that the two cases are markedly different for low neutrino 
luminosities (late times) . If all the initial average neutrino energies and luminosities were 
the same neutrino oscillations would clearly have no impact. Our analysis highlights 
the significance of having a precise knowledge of the muon and tau neutrino spectra and 
luminosities in assessing core-collapse supernovae as a site of r-process nucleosynthesis. 
In both of the two generic, but markedly different, cases we studied we found that 
neutrino self-interactions are crucial in setting the value of the electron-fraction to be 
rather large. We should emphasize that our treatment of the neutrino self-interactions 
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in Eqs. (I45|) and (j46p is an approximation using flux-averaged values. Hence an exact 
treatment of the neutrino self-interactions remains to be the key issue in the proper 
description of the neutrino transport in core-collapse supernovae. 
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